Search for: All records

Abstract. Peatlands have often been neglected in Earth system models (ESMs).Where they are included, they are usually represented via a separate, prescribed grid cell fraction that is given the physical characteristics of a peat (highly organic) soil. However, in reality soils vary on a spectrum between purely mineral soil (no organic material) and purely organicsoil, typically with an organic layer of variable thickness overlying mineral soil below. They are also dynamic, with organic layer thickness and its properties changing over time. Neither the spectrumof soil types nor their dynamic nature can be captured by current ESMs. Here we present a new version of an ESM land surface scheme (Joint UK Land Environment Simulator, JULES) where soil organic matter accumulation – and thus peatland formation, degradation and stability – is integratedin the vertically resolved soil carbon scheme. We also introduce the capacity to track soil carbon age as a function of depth in JULES and compare this to measured peat age–depth profiles. The new scheme is tested and evaluated at northern and temperate sites. This scheme simulates dynamic feedbacks between the soil organic material and its thermal and hydraulic characteristics. We show that draining the peatlands can lead to significant carbon loss, soil compaction and changes in peat properties. However, negative feedbacks can lead to the potential for peatlands to rewet themselves following drainage.These ecohydrological feedbacks can also lead to peatlands maintaining themselves in climates where peat formation would not otherwise initiate in the model, i.e. displaying some degree of resilience. The new model produces similar results to the original model for mineral soils and realistic profiles of soil organic carbon for peatlands.We evaluate the model against typical peat profiles based on 216 northern and temperate sites from a global dataset of peat cores.The root-mean-squared error (RMSE) in the soil carbon profile is reduced by 35 %–80 % in the best-performing JULES-Peat simulationscompared with the standard JULES configuration. The RMSE in these JULES-Peat simulations is 7.7–16.7 kg C m−3 depending on climate zone, which is considerably smaller than the soil carbon itself (around 30–60 kg C m−3). The RMSE at mineral soil sites is also reducedin JULES-Peat compared with the original JULES configuration (reduced by ∼ 30 %–50 %). Thus, JULES-Peat can be used as a complete scheme that simulates both organic and mineral soils. It does not requireany additional input data and introduces minimal additional variables to the model. This provides a new approach for improving the simulation of organic and peatland soils andassociated carbon-cycle feedbacks in ESMs. 

Abstract. In the global methane budget, the largest natural sourceis attributed to wetlands, which encompass all ecosystems composed ofwaterlogged or inundated ground, capable of methane production. Among them,northern peatlands that store large amounts of soil organic carbon have beenfunctioning, since the end of the last glaciation period, as long-termsources of methane (CH4) and are one of the most significant methanesources among wetlands. To reduce uncertainty of quantifying methane flux in theglobal methane budget, it is of significance to understand the underlyingprocesses for methane production and fluxes in northern peatlands. A methanemodel that features methane production and transport by plants, ebullitionprocess and diffusion in soil, oxidation to CO2, and CH4 fluxes tothe atmosphere has been embedded in the ORCHIDEE-PEAT land surface modelthat includes an explicit representation of northern peatlands.ORCHIDEE-PCH4 was calibrated and evaluated on 14 peatland sites distributedon both the Eurasian and American continents in the northern boreal andtemperate regions. Data assimilation approaches were employed to optimizedparameters at each site and at all sites simultaneously. Results show thatmethanogenesis is sensitive to temperature and substrate availability overthe top 75 cm of soil depth. Methane emissions estimated using single siteoptimization (SSO) of model parameters are underestimated by 9 g CH4 m−2 yr−1 on average (i.e., 50 % higher than the site average ofyearly methane emissions). While using the multi-site optimization (MSO),methane emissions are overestimated by 5 g CH4 m−2 yr−1 onaverage across all investigated sites (i.e., 37 % lower than the siteaverage of yearly methane emissions). 

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO₂ fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (V_cmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r² =  0.76; Nash–Sutcliffe modeling efficiency, MEF  =  0.76) and ecosystem respiration (ER, r² =  0.78, MEF  =  0.75), with lesser accuracy for latent heat fluxes (LE, r² =  0.42, MEF  =  0.14) and and net ecosystem CO₂ exchange (NEE, r² =  0.38, MEF  =  0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r² values (0.57–0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r² values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r² < 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized V_cmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average V_cmax value.